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ON THE TIME SLICING APPROXIMATION OF FEYNMAN 
PATH INTEGRALS FOR NON-SMOOTH POTENTIALS 

FABIO NICOLA 


Abstract. We consider the time slicing approximations of Feynman path in¬ 
tegrals, constructed via piecewice classical paths. A detailed study of the con¬ 
vergence in the norm operator topology, in the space B(L 2 (S . d )) of bounded 
operators on L 2 , and even in finer operator topologies, was carried on by D. 
Fujiwara in the case of smooth potentials with an at most quadratic growth. In 
the present paper we show that the result about the convergence in B(L 2 (Mr)) 
remains valid if the potential is only assumed to have second space derivatives 
in the Sobolev space H d+1 (R d ) (locally and uniformly), uniformly in time. The 
proof is non-perturbative in nature, but relies on a precise short time analysis of 
the Hamiltonian flow at this Sobolev regularity and on the continuity in L 2 of 
certain oscillatory integral operators with non-smooth phase and amplitude. 


1. Introduction 

Consider the Schrodinger equation 

(1) ihd t u = —~h 2 Au + V(t,x)u 

where 0 < h < 1 and the potential V(t, x), t e M, x £ W * l is a real-valued function. 
A fundamental problem is to find more or less explicit formulas for the solutions. 
As a general principle, R. Feynman in [9,10] conjectured that, under reasonable 
conditions on V(t,x), the propagator should be expressed as a certain formal in¬ 
tegral on an infinite dimensional space of paths in configuration space. While this 
principle was explored and made mathematically rigorous by several authors, here 
we will adopt the so-called “time slicing” approach as developed by D. Fujiwara 
and his school in [11-18,22,23,25-31,34,40,41], 

In short, inspired by the well-known formula of the Schrodinger propagator of 
the free particle, in general one considers the parametrices E^ 0 \t,s) defined by 

< 2 > s >'<*> = ,2 whw* L dv ■ 
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where S(t,s,x,y ) is the action computed along the classical path 7 (i.e. the path 
satisfying the Euler-Lagrange equations) satisfying the boundary condition 7 ( 5 ) = 
y, 7 (t) = x (which will be unique for 1 1 — s\ small, under the assumptions below). 
Hence 

(3) S(t, s, x,y) = J £( 7 (r), 7 (r),r)dr, 

where C is the Lagrangian of the corresponding classical system. 

Then one considers a subdivision U : s — t 0 < £ < ... < — t of the interval 

[s, t] and the composition 

(4) f? (0) (h, t, s ) = EW(t, t L _ 2 )... U (0) (U, s ), 

which should be regarded as an approximation of the path integral, exactly as 
the Riemann sums appear in the definition of the Riemann integral. Therefore, 
one wonders whether the operators E^ 0 \Q, t, s) converge to the actual propagator 
U ( t , s) as 

u;(fl) := sup{fj — tj-i : j — 1 ,..« t L} —> 0. 

The problem of the convergence of E^ (H,f, s) in the norm operator topology, in 
the space H(L 2 (R d )) of bounded operators in E 2 , was first considered in [ 11 ], where 
V (■ t , x) was assumed to have partial derivatives 

d*V G L°°(M x R d ) for |a| > 2 . 

Under the same hypotheses it was then proved in [12] that there is in fact con¬ 
vergence in a finer topology, at the level of the integral kernels of the involved 
operators. This suggests that the former result of convergence in H(L 2 (M d )) could 
hold under weaker regularity assumptions on V(t,x). The aim of this note is ex¬ 
actly to investigate this issue. 

We consider the so-called Kato-Sobolev spaces iL[) z (M d ), (also called uni¬ 

formly local Sobolev spaces), of functions / G L) oc (M d ) which belong to the usual 
L 2 -based Sobolev spaces H K (B) for every open ball B C R d of radius 1, uniformly 
with respect to B , in the sense that sup B ||/||h 7 B) < + 00 . In comparison with 
the usual Sobolev spaces H K (R d ), the functions in for k > d/2 are still 

bounded and continuous but need not decay at infinity; see [2,3,21] and Section 2 
below for details. 

Now, we assume the following condition. 

Assumption (A) V(t,x) is a real-valued function in L] oc f Rx M d ) such that for 
almost every t G M and |a| < 2 the derivatives dfV(t,x) exist and are continuous 
with respect to x, and moreover 

(5) dfV G L°°(M;i^ + 1 (M d )) for |a 


2. 
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Under this assumption we will see that there exists 5 > 0 such that if 0 < 
\t — s\ <5 and x, y G there exists only a classical path 7 satisfying y(s) = y 
and 7 (t) = x, so that S(t,s,x,y ) is well defined. Also, the operators E^(t,s), 
initially defined on Schwartz functions, extend to bounded operators on L 2 (W l ), so 
that the composition s) in (4) is well defined too, for any subdivision fl 

with uj(Q) < 6. 

We have the following result. 

Theorem 1.1. Suppose the condition in Assumption (A). For every T > 0 there 
exists C = C(T) > 0 such that for 0 < t — s < T and any subdivision of the 
interval [s,t] with ca(fi) < 5, and 0 < h < 1, we have 

( 6 ) ||£' ( 0 ) (n,f, s) - U(t,s)\\ L 2^ L 2 < Cu(Q)(t - s ). 

We notice that the rate of convergence of (Q,t: s ) is the same as that which 
appears in [ 11 , 12 ] for smooth potentials. 

Let us observe that usually non-smooth potentials, even more singular than 
those considered above, are successfully treated as a perturbation. For example, 
if V(x) is a time-independent potential in L°°(R d ) (or in dimension d = 3, V G 
L 2 (M 3 ) + L°°(M 3 )) one has the Trotter formula (cf. [36, Theorem X. 66 ], [38]) 

e"H tH = s- lim (e-£ Ho e -£ v Y 

n—>00 \ / 

where H = H 0 + V(x), H 0 = — (l/2)h 2 A. However, we should notice that 
the operators {er^ H °e~^ v \ n are in general different from the approximations 
Q,,t,s = 0) in (4) and, most importantly, the limit above holds only in the 
strong operator topology (the Trotter formula is also valid with norm convergence 
in some nontrivial circumstances, but for parabolic equations; cf. [19]). Instead the 
main point in Theorem 1.1 is the convergence in the norm operator topology. It 
seems that this issue has never been considered for unbounded and time-dependent 
non-smooth potentials. 

We also like to mention the approach to path integrals as infinite dimensional 
Fresnel integrals, developed by S. Albeverio and coworkers (see [1,33] and the 
references therein), which allows one to treat potentials of the type “polynomial + 
Fourier transform of a finite measure”. However, again, it does not seem that the 
authors consider the problem of the convergence in the norm operator topology. 

The choice of the Sobolev exponent d + 1 in (5) is justified as follows. All the 
known results about the continuity in L 2 (R d ) of oscillatory integral operators such 
as ( 2 ), with non-smooth phase (and amplitude) seem to require that the second 
derivatives of the phase have, roughly speaking, at least s additional derivatives 
both with respect to x and y, for some s > d/ 2 (when the derivatives are meant 
in the L 2 -Sobolev scale; see [3,37] and the references therein). This fact is at 
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least heuristically related to the known counterexamples about the continuity in 
L 2 of pseudodifferential operators (see e.g. [3, Remark 2]), since often the second 
derivatives of the phase play the same role in the estimates as the amplitude; cf. 
(40) below. Now, under the assumption (A) we will prove that S(t,s,x,y ) has in 
fact second space derivatives in H^~ 1 ( R 2d ), hence possessing d /2 + 1/2 additional 
derivatives both with respect to x and y. One could probably refine the result 
considering in (5) the fractional Sobolev space for some s > d, but we 

preferred to avoid further technicalities. 

In the proof we follow the same strategy as in the case of smooth potentials 
[11,12] but we need, as a preliminary result, a refined short time analysis of the 
Hamiltonian flow in such low regularity spaces. This represents, in fact, the main 
technical issue of the paper and is achieved via a priori estimates using, as tools, 
compositions and inverse mapping theorems for Sobolev mappings [4,20]. Also, 
the continuity results for oscillatory integral operators with non-smooth phase and 
amplitude proved in [3] will play a key role. 

We observe that, under additional space regularity for the potential, we can prove 
similar results for the convergence of higher order parametrices, where powers of h 
appear in the right-hand side of (6). This seems of particular interest in view of 
the semiclassical approximation and will be studied in Section 5. 

As a final remark let us mention a different yet non-perturbative approach, where 
the Schrodinger propagator is constructed by superposition of wave packets [7,8,24, 
32,39]. That approach would be certainly worth investigating in connection with 
path integrals, especially for classes of potentials which locally have the regularity of 
a function whose Fourier transform is in L 1 (M d ) (cf. [6,8]); however the parametrices 
in that case do not have anymore the oscillatory integral form in (2), and the 
problem of convergence would concern a different type of approximations. 

The paper is organized as follows. In Section 2 we prove some preliminary results 
about the composition and the inverse of Kato-Sobolev mappings, and recall a 
sufficient condition for the continuity in L 2 of oscillatory integral operators. Section 
3 is devoted to a detailed short time analysis of the Hamiltonian flow, the focus 
being on the Sobolev regularity issue. In Section 4 we prove Theorem 1.1. Finally 
in Section 5 we refine the above result for higher order parametrices, provided V 
itself is more regular. 


2. Preliminary results 

2.1. Sobolev spaces. We adopt the usual notation H K (B) = W K,2 (B) and H K (JR. d ), 
k6 N, for the L 2 -based Sobolev spaces on an open ball B C W l or on M d , respec¬ 
tively. We also consider the fractional Sobolev spaces FP(M d ), with sgR. 
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In the sequel we will often use the Gagliardo-Nirenberg-Sobolev inequality (see 
e.g. [35]) 

(7) IISYIGm < C||/Iiy^“ 

valid for \a\ < a, 1/p = \ol\/{2k), where B is any open ball of radius 1 in for a 
constant C > 0 independent of B. 

2.2. Kato-Sobolev spaces. For k G N we consider the space of functions 

/ in Lj oc (R d ) satisfying 

||/||h».(r<«) : = sup \\f\\H-(B) = sup sup \\d a f\\ L 2 {B) < +oo 
B B |q|<k 

where the supremum is made on all open balls B C M d of radius 1, and the 
derivatives are meant in the distribution sense. In fact, considering balls of any 
fixed radius r > 0 would give equivalent norms. 

More generally, for every s G R one defines the spaces H s ul {W l ) of temperate 
distributions / G S'(/R d ) such that 

\\f\\H° ul (Rd) ■= SUp ||x(- — y)f\\H*(R d ) < +00, 

y€ K d 

where \ ^ \ {0}. Changing y gives equivalent norms, and it is also easy 

to see that this definition reduces to that above when s is a non-negative integer. 
Moreover if s > d/ 2, C C(M d ) D L°°(M d ) is a Banach algebra [2,3,21], 

In Assumption (A) in Introduction we used mixed norm spaces which are defined 
precisely as follows: if n G N we denote by L°°(M; H/^/W 1 )) the space of functions 
f(t,x ) in L ] oc {M x R d ) such that their distribution derivatives d/f for |a| < k are 
in Lj oc (R x M . d ) and moreover 

II f \ \ L°° '■= ess-sup || f(t, •)||ij« i (Rd) < +CXD. 

We will need the following property about scaling. 

Proposition 2.1. Let kgN. With the notation a t (x,y ) = a(x,ty ) for 0 < t < 1, 
x, y G R d , there exists a constant C > 0 independent oft and a such that 

ll a i|lH« i (R 2d ) < C'||a||^ + [d/2]+i^ ]K2d ^. 

Here [d/2] stands for the integer part of d/2. 

Proof. Let |a| + |/3| < k and B be an open ball of radius 1 in M. d . We have 

\\dfd^a t (x,-)\\ L 2 {B) = t m - d/2 \\dfd^a(x,-)\\ L 2 {§) 

where B = {ty : y G B} is a ball of radius t. Now, if |^| > d/2 we have t^~ d / 2 < 1 
and one concludes easily by taking the norm in L?{B' ) of this expression, where 
B' C is any ball of radius 1. 
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If instead \/3\ < d/2, hence \/3\ < m : = [d/2] + 1 we set 1/p = |/?|/(2 m) and we 
continue the above estimate, by Holder’s inequality, as 

< ct'“ n i -'</< 2 ”>» n^afo(i, •)||l, ( b») 


where B" D B is the ball of radius 1 with the same center as B. By the Gagliardo- 
Nirenberg-Sobolev inequality (7) this last expression is dominated by 


I H m (B") 


which in turn is dominated (since m > d/2, and therefore d /( 2m )) < 1 and 
H m (B") -A L°°(B")) by 

C"\\d%a(x, -)|| h™(b”), 

and one concludes as above. ■ 


Remark 2.2. One easily sees that a loss of d/2 derivatives in Proposition 2.1 is 
inevitable. In fact, suppose that, for \a\ = n and some r > 0 the following estimate 
holds: 

\\dfat\\ L 2(BxB) < C'||a||^+q K 2d) 

where B = {x G R d : |a;| < 1}. We test this estimate with a{x,y) = aW (\x)a/ 2 \\y), 
\ = 1/t >1, where af 2 \y) are in Cf°(B) \ {0}. We obtain the estimate 

||<9 a (a^(A-)) || l 2 (b) ||° (2 ' ) ||l 2 (b) < G'lla^^A^a^^A^II^/c+r-^d) 

and therefore as A —>■ +oo we deduce \ K ~ d / 2 < C"\ K+r ~ d , which implies r > d/2. 

We are particularly interested in the issue of the composition and inverse map¬ 
ping theorem for Sobolev mappings. 

Proposition 2.3. Let g :R d —>■ M cZ be a globally bi-Lipschitz map, i.e. satisfying 
(8) Cf^x-y] < \g(x) - g(y)\ < C 0 \x - y\, x,y eR d 

for some constant Co > 0. 

Let n G N. Then, if f e D L°°(R d ) and Dg e R dxd ) we have 

f o g G D L°°(R d ). 

Here Dg stands for the Jacobian matrix of g. 

Proof. The result is a variant of the known composition formulas for Sobolev map¬ 
pings (see e.g. [4,20] and the references therein). We provide a short proof for the 
benefit of the reader, since we did not find the statement exactly in the present 
form in the literature. 
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It is clear that for every open ball B of radius 1, 

11/ o g\\mB) <C\\fo g\\ L oo( B ) <C\\fo 5 f|| L oc {K d) < C'||/|| L oo (K d ) . 

To estimate the derivatives of / o g we use the chain rule (for the justification of 
the chain rule at this regularity we refer to [4]). Namely, let g — (<y 1? .,. ,gd)- We 
can write d a (f o g) ) 1 < |a| < k, as a linear combination of terms of the form 

d a f(g(x)) d^gj^x)... d^g jH (x) 

where 1 < |cr| < |a|, \gi \,..., \g\ a \\ > 1, pi + ... + g\ a \ = a, and j\,... ,j\ a \ G 

{l,...,d}. 

Now, given an open ball B C of radius 1, we estimate the norm in L 2 (B) of 
this expression by Holder’s inequality as 


( 9 ) II (9 a f) o g\\LPo(B)\\d^ gn\\LPi(B) ■ ■ ■ Wd^gj^ \\ l p w\(b)i 

where we choose 


o' 


|/ij 


Po 2 K ’ Pj 


Observe that, in fact, 


o r 


ki 

yi 

% P, 2 K 


2k 


a - cr 


3 = 


2k 


Kb 

2k 


< 


On the other hand, by the bi-Lipschitz assumption (8) we have, with B = g(B), 

\Wf)og\\ LV0{B) <CWf\\ LP0( By 

Since B = g(B) can be covered by a fixed number of balls B' of radius 1, and by 
the Gagliardo-Nirenberg-Sobolev inequality (7) we can continue the estimate as 

II (d a f) o g\\ L p 0 ( B ) < C' sup \\d a f\\ LP0 ( B >) 

B' 


< C" sup 

B' 


l-\a\/K,\\ r ||kl/K 
L°°(B')\\J II H K (B') 


< ci/iik^ii/iiSk 

Hence we obtain 

sup || (d a f) og \\ L po( B ) < +oo. 

B 

Similarly we estimate the other factors in (9): for j = 1,..., |oj we have | fXj\ > 1, 
so that by the Gagliardo-Nirenberg-Sobolev inequality we obtain 

which is uniformly bounded with respect to B by the assumptions on g (observe 
that (8) implies Dg G L°°). ■ 
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We also recall the following inverse mapping theorem. 

Proposition 2.4. Let seN. Let g : W l —> be a globally bi-Lipschitz map, i.e. 

satisfying (8), with Dg G Hff l {U. d ]M dxd ). Then Dg- 1 G i/^(M d ;R dxd ) as well. 

Proof. By applying the inverse mapping theorem for Sobolev mappings as stated 
e.g. in [4, Theorem 1.1] (with m = /c+1 and p — 2) one obtains that Dg- 1 G H K (B) 
for every open ball B C even if g were only locally bi-Lipschitz. Actually, since 
g is globally bi-Lipschitz and Dg is in H K locally uniformly, an inspection of the 
proof of [4, Theorem 1.1] shows that the desired estimates are uniform with respect 
to the balls B, provided these latter have the same fixed radius. This concludes 
the proof. 


Remark 2.5. It also follows from the proof of Proposition 2-4 that we have, in fact, 
\\Dg^ 1 \\H K < C for a constant C depending only on Co, k, d and upper bounds for 
\\Dg\\n :i , where C 0 in the constant in (8). A similar remark applies to Proposition 
2.3. In the sequel we will apply freely these “uniform versions” of Propositions 2.3 
and 2.4 without further comments. 


2.3. Oscillatory integral operators with non-smooth phase and ampli¬ 
tude. Let A be an oscillatory integral operator of the form 

(10) Af(x) — f e lXS{x ' v) a(x,y)f(y)dy, A > A 0 > 0. 

JR d 


Several conditions on the phase function S(x,y ) and the amplitude a(x,y ) are 
known for A to be bounded on L 2 (M d ); see e.g. [37] and the references therein. 
Here we recall the following result from [3, Corollary 1], 


Proposition 2.6. Let A be the operator in (10), and suppose S real-valued with 
d a S G H^ ll (R 2d ) for |a| = 2, and a G H* l (Ml 2d ) for some s > d. Assume, moreover, 
that 


( 11 ) 


, d 2 S . , 

det—^— {x,y) 
oxoy 


> 5 


x, y G R d 


for some 5 > 0. 

Then A, initially defined on Schwartz functions, extends to a bounded operator 
on L 2 (R d ). In fact there exists a constant C > 0 independent of A > Ao, S, S and 
a such that 

PllL»-y < Cr I A-« 2 exp(C||D 2 5|| H . I )||a|| H . I , 


Here D 2 S denotes the Hessian matrix of S. 
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3. Short time analysis of the Hamiltonian flow 

3.1. Sobolev regularity of the Hamiltonian flow. In this section we will as¬ 
sume the following hypothesis on the potential V(t,x). 

Let k G N, k > d + 1. 

Assumption (B) V(t,x ) is a real-valued function in L) 0C (IR. x M d ) such that for 
almost every tel and \a\ < 2 the derivatives dfV(t,x ) exist and are continuous 
with respect to x, and moreover 

dfV G L°°(M; (M d )) for |a| = 2. 


In particular, we will take k — d + 1 in Theorem 1.1 and higher values of k in 
Theorem 5.2 below. 

Denote by (x(t, s, y, r]),f{t, s, y, y)), s,t G M, y, rj G the solution of the Hamil¬ 
tonian system 

(12) x = €, £ = ~V x V(t,x) 

with initial condition at time t — s given by x(s, s, y, rj) = y , £(s, s, y , rj) = rj. 

Since HjjW 1 ) C L°°(R d ) (because k > d + 1 > d/2), for |a| = 2 we have 
dfV G L°°(M x W l ), and moreover dfV(t,x) is continuous with respect to x for 
almost every t, so that the flow 

(13) (y, rj) eA (x(t, s, y, rj),£(t, s, y, ij)) 

defines a C l canonical transformation W l x M, d —y M, d x W l . Moreover, for fixed s, as 
a function of t, y, rj, the function x(t, s, y, rj) is C 1 whereas the function j[t, s, y, rj) 
is locally Lipschitz. 

Concerning quantitative information, we observe that it follows at once from 
Gronwall’s inequality that for any T > 0 there exists a constant C = C(T) > 0 
such that 


(14) 


dx 

dy 


(t,s) 


+ 


dx 
d rj 


M 


d£(. \ 


+ 


d£(. ) 


<C, 


if \t — s| < T. In particular the map (x(t,s),£(t,s)) is globally Lipschitz. 

In fact, we can prove the following result, which is finer when \t — s\ is small. 

Lemma 3.1. For every T > 0 there exists a constant C = C(T ) > 0 such that, 
for \t — s\ < T, 


(15) 


dx 

dy 


(' t , s) — Id 


< C\t — s| 2 , 


dx 

drj 


(t,s) 


< C\t — s\ 
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(16) 

oy 11 z, 00 (Kg xu,,. 

Proof. From (12) we have 


<C\t-s\, ||||(t, s) - Id 


L°°(R^xK^) 


< C\t — s| 2 . 


^ {t ' s) = ^ + lpf s)dT 


where u stands for yj or pj, j = 1,..., d. Hence we obtain 

A, 

du 

rt d 2 


% (t ' a) -Tu= I ^ (T ' s)dT 


and 


df / \ dy 

SI (M) “aI + 


axdx VlT - x<T ' s)) ^ dT 


t pr 



P^-V(r,x(r,s)) d f- ( m )*. 


Then, using (14) we obtain at once the desired estimates. 


Our aim now is to prove that the flow (x(t, s, y, y),f(t, s, y, y)) has the same space 
regularity as V. X H, namely that the first space derivatives of x(t, s, y, y),f(t, s, y , y) 
belong to Hf Ll (Mi 2d ). In fact, the following precise asymptotic estimates of their 
Sobolev norm as \t — s\ —y 0 will play a key role. 

Proposition 3.2. Let T > 0. There exists a constant C > 0 depending only 
on T and upper bounds for the norm of D 2 V in L°°(M; such that, for 

2 < |a| + \/3\ < k + 1 ? 

(17) S )|U= (B) < C\t - S |'"' +2 , ||^(«, S )|| 12 ( B ) < C\t - 

for every open ball B C M, d x M. d of radius 1. 

Here we used the notation D 2 V for the Hessian of V with respect to x. 

Proof. First Step: Regularization. 

First of all we observe that it is sufficient to prove the desired estimates under 
the additional assumption that V ( t , x) is smooth with respect to x for almost every 
tel. The argument uses standard techniques, but we provide a sketch for the 
benefit of the reader. 

Assume that the result holds for smooth potentials, and consider smooth regu¬ 
larizations 

v e (t, ■) = V(t, ■) *pe, 0 < e < 1 , 
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where p e (x) = e d p(e 1 t) is a standard modifier in M d , i.e. p G C^°(W l ), p > 0, 

h*P( x ) = L 

We have || D 2 V e (t, -)||< C\\D%V(t, •) 11 h* :1 for a constant C independent of e, so 
that the corresponding solution (x e (t,s),£ e (t,s)) will satisfy the estimates in (17) 
with a constant C independent of e. 

On the other hand, we easily verify, in this order, the following facts. 

To simplify notation, set X = (x, £), b(t, X) = (£,,—V X V(£, x)), X e = (x e ,£ e ), 
b e(t,X) = (£, -V x V e (t,x)),Y = (y, V ). 

a) V x y(. s 0) G Lj oc (M.) and there exists a constant C > 0 such that 

\X x V e (t,x)-X x V(t,0)\ <C(l + |x|) 

for almost every t G R and every x G R d , 0 < e < 1. 

Indeed, for almost every iGlwe have the Taylor expansion 

d 

(18) V(t,x) — a Q (t) + J2 a i(t) x j + v(2 \t,x), 

3 =1 

where I/ (2) G L“ C (R x R d ), because D 2 I7 G L°°(R x R d ) and D 2 V{t,x) is 
continuous with respect to x for almost every t. Since V G Lj oc (R x R d ), we 
see that the function ao(t) + 'Yhj=\ a j{t) x j is in Lj oc (R) for almost every x, 
which implies aj G Lj oc (M), j = 0,..., d. In particular V X V (•, 0) G Lj oc (R). 
Similarly we have 

V X V (t, x) - X x V(t, 0) = b(£, x) 

where |b(t, x)| < C |x| for almost every 1 G R and every x G R d . Since 
V x y e (£, •) - V X F (t, 0) = b(t, •) * p e 
the desired conclusion follows easily. 

b) The solutions X e (t, s, Y ) for fixed s G R are bounded on the compact subsets 
ofR x R 2d , uniformly with respect to e. 

In fact, since X e (t,s,Y) = (x e (t, s,Y), £ e (t, s, Y)) is the flow of b e (£, X) 
we have 

X e (t,s,Y) = Y+ j (0, —V X V (t, 0))dr 

J S 

+ [ (£ £ (t, s, Y), —V x K(r, x £ (r, s, F)) + V X V(t, 0))dr. 

J S 

The estimates in the previous point and Gronwall’s inequality give the de¬ 
sired conclusion. 

c) For almost every f G 1, b e (£, X) converges to b(£, X) as e —* 0, uniformly 
on the compact subsets of R 2d . 
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d) The difference b e (t, X) — b(f,X) is essentially bounded on the compact sub¬ 
sets ofM. x M. 2d , uniformly with respect to e. 

In fact with the notation in the point a) we have b G L“ C (K x R d ). On 
the other hand 

v,K(t, •) - V x V(t, •) = b (t, -)*p e - b (t, •) 

and the claim follows at once. 

e) For fixed s6l, X e (t,s,Y) converges to X (t,s,Y) as e —» 0, uniformly on 
the compact subsets o/RxR d . 

In fact it turns out that for fixed T,R > 0, by the point b) there exists 
a ball B' C M? d where the functions X1 e (t,s,Y) take values for |f| < T, 
\Y\ < R, and we have 1 

\X e (t,s,Y) - X(t,s,y)| < ||b e (r, •) — b(r, •)||l°°(b') dr 

+ j % || Vxb(r, 0 IIloo(r M) | X e (r, s, Y) - X(r, s, Y)\ dr. 

The first integral tends to zero by the dominated convergence theorem (by 
the points c) and d)) and one concludes by Gronwall’s inequality. 

Hence for every open ball B C R 2d the functions X e (t,s,-) = (x e (t,s),£ e (t,s)) 
converge in V\B) to X(t,s) = (x(t,s),£(t,s)) as e —* 0, and then the desired 
estimates (17) hold 2 for (x(t,s),£(t,s)). 

Second step: A priori estimates. 

We now prove the formulas (17) as priori estimates. In order to apply an in¬ 
ductive argument we have to formulate the inductive hypothesis in the following 
stronger form: for 

2 < |a| + \/3\ < k + 1, \t - s\ < T, 1/p = (|a| + \j3\ - l)/(2/c), 


we have 


(19) 

\\ap^x(t,s)\\ Lrm <c\t-s\^ 

(20) 

|| 9 “aJqt,s)l| W (B) < c\t - s| |fi|+1 


for a constant C > 0 independent of B, as in the statement. Therefore we prove 
these estimates by induction on |a| + \f3\. 

1 The L°° norm of a vector-valued function / is meant as the L°° norm of the Euclidean norm 
1 / 0 * 01 - 

2 As a conseciuence of the Riesz representation theorem and the density of Cf°(B) in L 2 (B), if 
u n is a bounded sequence in L 2 (B), say |jit„||x,2(s) < C, converging to a distribution u in T>\B), 
it turns out u £ L 2 (B) and ||u||l 2 (b) < C. 
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We have, for 2 < \a\ + \(3\ < k + 1, 


( 21 ) 

and 


dyd^x(t, s, y, rf) = 



dydfe{r,s,y,'ri) dr 


(22) d*dfe(t,s,y, r)) = - I ^{d Xj V x V)(t, x(t, s))d^d^Xj(T, s, y, rj) dr 

ds j = 1 

- [ b a A T ’S,y,v)dT 


where b a A r i s ,y,v) is a linear combination of terms of the form 



where j u ... ,j\ a \ e {1d}, v x + ... + v\„\ = a, H\ + ... + n\ a \ = (3, \vj\ + \pj\ > 1 
for j — 1,.,. , \cr\ and 2 < \cr\ < |a| + \/3\. 

We can estimate the norm in L P (B) of this expression, with 


by Holder’s inequality as 


1 |a| + \/3\ — 1 

p 2 k 


(23) ||(t?V 1 ,l/)(T,i(r, S ))|l t ,„ (B) ||S^S«i|| 1 „ (B) ... \\dy'°%'"'x 


\L p \°\(B)i 


where 

1 

Po 

Observe that 


I - 1 1 = M + \di\ ~ 1 

2k ’ pj 2 k 


H 1 

E- 


a | — 1 + |ck| + \f3\ — \a\ 1 

2k p 


Since the map (x(t,s),£(t,s)) is a canonical, therefore measure preserving trans¬ 
formation, the first factor in (23) can be written as 


||(^VW)(r,i(T,i,))||„ (B| = ||(<?VW)(t,.)|| 1 „o ( b) 

where B C x is the image of B by the flow (x(t, s), £(t, s)) (we are thinking 
of V(t,x ) as a function of r,x,£, constant with respect to £). Now, using (14), 
we see that B can be covered by IV = N(T,d ) boxes B’ x B”, with B',B" C 
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balls of radius 1, for \t — s| < T, so that we can continue the estimate, by the 
Gagliardo-Nirenberg-Sobolev inequality (7) (using 2 < \a\ < n T 1) as 

\\(d"\7 x V)(T,x(T,s))\\ L r 0{B) < Csup\\d"V x V(T,-)\\ L P 0{B ») 

B" 

< r'snn II D 2 I/ll 1_ 6°'T 1 )/'q| n 2 T/||(l°'l^ 1 )/ K 

— b ^P\\ 1 J x v IIL°°(RxK d ) W 1 J x v IIl°°(R;IJ''(B")) 

/ (~ii \I 7-)2t / ||1-(|ct|-1)/k|| 7-)2t am (|ct|—1)/k 

— \\- Ly x v llL°°(RxR d ) \\ 1 J x v llL°°(R;ff« i )- 

The last expression is finite, since 

II Ar^ll-L°°(RxR d ) < CWDIVWloo^.r^) < TOO 

by assumption. 

To treat the other factors in (23) observe that 

\\ar‘d»x(T, s)|| m , B) < C|r - s|l«i 


for a constant C = C(T ) if |r — s\ < T. This holds for \uj\ + \nj\ = 1 by (15) and 
if \uj\ T \fJ>j\ >2 by the inductive hypothesis (because \uj \ + \fij\ < \a\ T \/3\ — 1, 
therefore in the latter case we have \a\ T \(3\ > 3). 

In conclusion, since Hi + . . . T /x \„\ = (3 we have 


\\ba,p(T, s)\\ LP{B ) < C\t - s | l/31 
and we can estimate in (22) 



This estimate, together with (21) and Gronwall’s inequality gives (20) and also 
(19). ■ 


Consider now the map 

(y, 0 I-T (x(t, s, y, C s, y, 0 ) := (x(t, s, y, (/{t - s)), (t - s)£(t, s, y, (/(t - s)), 

with |/,(e W l . Observe that it is a C 1 canonical transformation. 

We have the following result. 

Proposition 3.3. We have, for j, k — 1,..., d, 

dec ■ 

(24) -Q- = {t s) a jk (t, s, y, £), 

(25) = s j,k~ (t- s)' 2 b jk (t,s,y, C), 

dx • 

(26) -y- = 8 jjk - (t - s) 2 c jk (t , 5, y , c), 
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(27) 1^- = (t- s) 2 c' jk (t,s,yX), 

where aj k (t, s ), bjk(t, s ), Cj k (t , s ), c' fc (£, s) are functions in a bounded subset o/i/^(R 2d ), 
if 0 < \t — s\ <T, for every fixed T > 0. 


Proof. Let us prove formula (25); the other formulas can be deduced similarly. 
First of all we observe that (25) defines the functions bj k (t,s), since Now, if 


B C 


ly x R^ is any open ball of radius 1 we have clearly by the second formula in 
(16) that, with r/ = (/(t — s), 


\\b jk (t,s)\\ L 2 {B) — (t s 
< C(t — s 


-2 ii d£j 


"dr) k 

-2 ii 


dr/k 


s)) - S jk \\ L 2 {b) 

(t, s, y, C/(t ~ s)) - Sjk\\L°°(B) < C 


for 0 < \t — s\ < T and some C' = C'(T) > 0 independent of B. 

We now estimate the derivatives of bj k . For l<|a| + |/3|<Awe have 

dyd(b jk (t, s, y,C) = ~(t~ (d%d% |^) (t, s, y, (/(t - s)) 

and therefore 

where 

B = {(y, V )e R d xR d : (y, {t — s)r]) G B}. 

Observe that B can be covered by N = 0(\t — s| _d ) balls B' C R d x R^ of radius 
1 , and therefore 


<c\t-s\ d/2 swp\\d%dP^(t,s)\\ L 2 {BI) < C'\t — s\ rf/2+l/3|+2 , 


d£j 


where in the last inequality we used the second formula in (17). 

Summing up we have 

\\dPdlb jk {t,s)\\ L 2 {B) < C 

for some C = C(T ) > 0 independent of B , for 0 < \t — s\ < T. ■ 


It follows from the proof, or more simply from the inclusion R 2d ) C L°°(M. 2d ) 
(recall k > d + 1), that the functions aj k (t, s),bj k (t, s ), Cj k (t, s), c' fc (f, s) belong to a 
bounded subset of L°°(R 2d ) for 0 < \t — s\ < T. 

In particular, we have obtained the following regularity result. 
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Proposition 3.4. For any fixed T > 0, the components of the canonical transfor¬ 
mation (x(t,s,y,(),f(t,s,y,()), have first space derivatives in a bounded subset of 
when 0 < \t — s\ < T. 

The now study the regularity of the inverse function of ( H» x(t, s, y, (), for \t — s\ 
small. 


Proposition 3.5. There exists S > 0 such that for 0 < \t — s\ < <5 and y G R d the 
map 

C ha x(t, s, y, C) = x(t, s, y, (/(t - s)) 

is invertible W 1 —» R d mid its inverse ( = ((t,s,x,y) has first derivatives with 
respect to x and y in a bounded subset of Hf d (R 2,i ). 

More precisely, for j, k — 1,.,., d, we have 
fjr. 

(28) S ’ V ^ = 5jk s ’ y) 

and 

(29) |^(C s, x, y) = -5 jk - ft - sfd' jk (t , s, x, y) 

where djkit, s,x,y) and d'- k {t, s,x,y) belong to a bounded subset of Hffi Rf x R]j) 
for 0 < \t — s\ < 5. 

Proof. Since Hfffi R 2rf ) is a Banach algebra (a > d + 1), by (24) we have 

dx 

(30) det —(t, s, y, () = 1 - (t - s) 2 a(t, s, y, () 

for some function a(t,s,y,Q in a bounded subset of Hffi R 2d ) C L°°(R 2d ), for 
0 < \t — s| < T, for any T > 0. Choose 5 > 0 such that 

(31) (t — s) 2 ||a(f, s, ■, -)||l°° < 1/2 


for 0 < \t — s\ < S. 

The invertibility of the C 1 map Q ha x(t, s, y, (), for 0 < \t — s\ < 5, then follows 
from Hadamard’s global inversion theorem [5, Theorem 2, page 93]. Moreover, 
for 0 < \t — s\ < 5 the map (y, () ha (x(t, s, y, (), y) will be globally bi-Lipschitz, 
uniformly with respect to t, s, with first derivatives in a bounded subset of Hfffi R 2a! ) 
(by (24)). By Proposition 2.4 we see that the same holds for the inverse map 
(x,y) ha (y,C(t,s,x,y)). 

In order to prove (28) and (29) we use the formulas 


d( 

dx 


(t, s, x, y) 


dx 


with C = C (t,s,x,y) 
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and 

0(^ Ox _ — 1 ()x 

— (t, s, x, y) = - j — (t, s, y, C)J s, y, (), with C = ({t, s, x, y ). 

Now we claim that it is sufficient to prove for the inverse matrix the formula 

( 32 ) IdC^’ Sl Vl = 5jk ~ ^ ~ s ’ y ' Okfc=i,-,d 

with d" k (t, s, y, () belonging to a bounded subset of H^(R 2d ). In fact, since H^(R 2d ) 
is a Banach algebra, by (32), (26) and the composition property in Proposition 2.3 
(with the map (x,y) ha (y,C(t, s,x,y)) playing the role of the map g) we see that 
(28) and (29) follow. 

To prove (32) we observe that, by (30) and (31) we have 

( det -(M,i/,0) = 1 + (t - s) 2 a'(t,s,y,(), 


where 


a'(t,s,y, C) 


a(t,s,y, C) 

1 - (t - s) 2 a(t,s,y, C) 


is easily seen to belong to a bounded subset of if^(K 2d ) for 0 < \t — .s| <5 (one 
can use the chain rule and arguments similar to those in the proof of Proposition 
2.3, using k > d + 1). This together with (24) gives (32) and concludes the proof. 


3.2. Sobolev regularity of the classical action. We close this section with 
an analysis of the regularity of the action S(t,s,x,y ) defined in (3) for, say, 0 < 
\t — s\ < 6 where 6 > 0 is the constant appearing in Proposition 3.5. 

As already observed, under Assumption (B) at the beginning of this section, for 
| a | = 2 we have d°V G L°°(R x R d ) and d%V(t,x) is continuous with respect to x 
for almost every tel, which implies by the arguments in [11, Section 2] that the 
for every fixed s, the function S(t, s, x, y) is almost everywhere totally differentiable 
with respect to (t, x, y) and satisfies the Hamilton-Jacobi equation 

d t S(t, s, x, y) + ^\V x S(t,s,x,y)\ 2 + V(t,x) = 0. 

Moreover S(t,s,x,y ) is of class C 2 with respect to x,y and satisfies the equations 

dS 

(33) — (t, s, x , y) = £j(t, s, y, rj(t , s, x, y)) 


and 

(34) 


— it, s, x, y) = s, x, y) 
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lor - 1 ..- d. 

Proposition 3.6. For 0 < \t — s\ < 5 we have 

(35) S(t, s , x , j/) = + (t - s)u(t , s, x, y)$ 

where the functions dfd^u>{t,s,x,y), for |a| + \/3\ = 2, belong to a bounded subset 
ofH^R™). 

Proof. Consider for example the derivatives of S with respect to x. Using the 
formula 

u(t, s, x , y) = (t - sy 2 ( t - s)S(t , s, x,y)--\x- y | 2 

together with (33), (34), (25) and (28) we deduce easily (cf. [11, Formula (2.12)]) 
that 


dxndxk 


u{t,s,x,y ) = -b jk (t,s,y,((t,s,x,y)) - d jk (t,s,x,y) 


+ (t - s) 2 E bjm (t, s, y,c(t, s, X , y))dmk(t, s, x, y). 


where the functions bj k , dj k are defined in (25) and (28). This expression belongs 
to a bounded subset of Hf l (lBL 2d ) 1 for 0 < \t — s\ < 5, by Propositions 3.3, 3.5 and 
Proposition 2.3 (with the map (x, y) ha (y, ({t, s, x, y)) playing the role of the map 

9 )• 

Similarly one can treat the other second derivatives of S, using the formulas 


dxjdy k 


>(t,s,x,y) = (t - s) 2 - ^(t,s,x,y) + 6 jk = d jk (t, s,x,y) 


dyjdy k 


r d(- i 

’(t,s,x,y) = (t- s)- 2 [ - -^(t,s,x,y) - S jk j = d' jk (t, s,x,y) 


where the functions dj k (t, s,x,y) are defined in (29). 


4. Analysis of the parametrices and proof of Theorem 1.1 

In this section we suppose that V ( t , s) satisfies the hypothesis in Assumption 
(A) in Introduction, which corresponds to Assumption (B) in the previous section 
with k — d + 1. In particular all the machinery of the previous section applies with 
this value of k. 

First of all we see that, by Proposition 3.6, the operator E^°\t,s) defined in (2) 
satisfies the assumptions of Proposition 2.6, with s — d + 1 , A = h~ x /\t — s|, for 
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0 < \t — s\ < 5, possibly for a smaller value of S. Hence it extends to a bounded 
operator in L 2 (R d ) and verifies 

(36) \\E(°\t,s)\\ L ^ L2 <C 
for 0 < \t — s\ < 5. 

We now show that s ) converges strongly to the identity operator as t —> s. 

Proposition 4.1. For every f G L 2 (M a! ) we have 

lim E {0) (t,s)f = f 

t^s 

in L 2 (R d ). 

Proof. First of all we observe that, by (36), it is sufficient to consider / G S(R d ). 

The result is clearly related to the stationary phase principle, but it is not easy 
to justify its application at this Sobolev regularity. Instead we argue “by density”, 
since we already know from [11, Proposition 4.3] that the result holds when V 
satisfies 

(37) dfV G L°°{R x R d ) for |a| > 2. 

Now, let V f be smooth regularizations of V, as in the proof of Proposition 3.2. 
It is easy to see (using i/^ +1 (M d ) C L°°(R d ) and L 1 (M d ) *L°°(R d ) C L°°(R d )) that, 
if V satisfies Assumption (A) in introduction, the potentials V f enjoy the property 
in (37). 

Let S e (t,s,x,y ) be the corresponding generating functions and 

< 38) : = (2,n( t -. W ^ Sj K ~' S '^ y)f(V)dV - 

As already observed in the proof of Proposition 3.2, Df.V e (t, •) belongs to a bounded 
subset of H d ^ 1 (R d ), so that all the results in Section 3 hold (with n — d + 1) 
uniformly with respect to e. In particular, the functions S e (t, s,x,y) are defined 
for 0 < \t — s\ < 6 for some fixed <5 > 0, independent of e. 

Let us now prove that Ee°\t, s)f(x) converges to s)f(x) for every x G R d , 

if 0 < \t — s| < 5. 

We first claim that 

S e (t, s, x, y) —y S(t, s, x, y) pointwise as e —> 0 
if 0 < \t — s\ < 6. In fact we have, with obvious notation, 

f l i 9 

S e {t,s,x,y) = J —\Ce(r,s,y, y e (t, s,x,y))\ 2 - V e (T,x e (T,s,y,r} e (t,s,x,y)))dT. 

We then apply the dominated convergence theorem: to check the convergence 
of the integrand function for almost every r G R one uses the point e) in the 
proof of Proposition 3.2, the fact that r} e (t,s,x,y) —>• rj(t,s,x,y) pointwise (which 
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follows easily by contradiction) and the fact that V e (r,-) V (r, •) uniformly on the 
compact subsets of W l , for almost every r. To check that the integrand function 
is conveniently dominated for fixed t, s, x, y , it is useful to consider the Taylor 
expansion (cf. (18)) 

d 

V(t,x) = a 0 (t) + J2aAt) Xj + V^(t,x), 

3=1 

where a,j G Ll oc (R d ), j — 0,..., d, and V^> G L^ C (M x R d ). If the modifier p(x) is 
e.g. an even function when restricted to every coordinate axis, we have 

d 

V e (T,x) = V(t,x) * pe{x) = a 0 (r) J r S ^a j (r)x j + V {2) (t,x) * p e (x). 

3=1 

Hence the desired claim follows easily from the point b) in the proof of Proposition 
3.2. 

Now, by the claim just proved we have indeed 

E f\^s)f{x) -A E {0 \t,s)f(x) 

for every x G M d , as e —> 0, by the dominated convergence theorem, because 
f eS(R d ) C L\R d ). 

By the Fatou theorem and the definition of “lim inf”, for every p > 0 there exists 
e 0 > 0 such that 

II E w (t,s)f - f || L 2 < lim inf \\E^ 0) (t, s)f - f\\& 

e —>0 

< l|£ ; i 0 0 | («.s)/-/llw+A«. 

Now, we know from the analogous result for smooth potentials in [11, Proposition 
4.3] that 

=o 

and therefore we conclude that 

lim sup ||TP (0 - > (t, s)f — /||z ,2 < p. 
t—ys 

Since p is arbitrary, we obtain the desired conclusion. ■ 

We continue the study of E ^ (t, s ) by observing that it is a parametrix in the 
sense that 

(ihd, + lft 2 A -V(t,x))E ( 0 \t,s)f = G m (t,s)f 


( 39 ) 
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with 

(40) 


G (0) (*,»)/ 


(l/2)ih(t — s) 
(2ni(t — s)h ) d / 2 


A^(t,s,x,y) f (y)dy, 


cf. [12, Formula (1.12)], where u>(t,s,x,y ) is defined in (35) (differentiation under 
the integral sign is justified e.g. for / G <S(R d ) exactly as in [11, Proposition 4.5]). 

It follows again from Propositions 3.6 (with k — d+ 1) and Proposition 2.6 (with 
s — d + 1) that G°(t, s ) is bounded on L 2 (R d ) and satisfies the key estimate 


(41) 


\\GW(t,s)\\ L ^ L2 <Ch\t-s 


for 0 < \t — s\ < 5, possibly for a smaller value of S. 

Concerning the actual propagator, we can easily prove its existence for more 
general potentials. 


Proposition 4.2. Let V(t, x ) be a real-valued function in Lj oc (R x R d ) such that for 
almost every lei and |a| < 2 the derivatives dfV(t,x) exist and are continuous 
with respect to x, with dfV G L°°(M x R d ) for |a| = 2. 

Let sGl. Then the Cauchy problem 

i ihd t u = —^h 2 Au + V ( t , x)u 
|u(s,t) = u 0 (x) 

is forward and backward globally wellposed in L 2 (R d ). 

If we denote by U (t, s ) the corresponding propagator, for every T > 0 there exists 
a constant C = C(T ) > 0 independent of h (0 < h < 1), such that 

(42) \\U(t,s)\\ L 2 ^<C 
for |t — s\ < T. 

Proof. Let y(£) be a smooth function in R d supported in |£| < 2, with y(^) = 1 for 
|£| < 1. We split the potential as 

(43) V{t,x) = X {Dx)V{t,x) + { 1 - X (Dx))V(t,x) . 

Vo (t,x) V\ (t,x) 

We claim that 


(44) dfV 0 G L°°(M x R d ) for |a| > 2 

and 


(45) Vj G L°°(R x R d ). 

Let us prove (44). For |a| > 2, write a = f3 + 7 , with |y| = 2 and 

dfV 0 (t,x) = x(D x )d%(dfV(t,x)). 
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We have dfV G L°°(M x M d ) by assumption and the Fourier multiplier x(D x )d% 
is bounded in L°°(R d ), because its symbol is the smooth compactly supported 
function 

Concerning (45), we write 

V 1 (t,x) = (l-x(D x ))A~ 1 A x V(t,x). 

By assumption A X V G L°°(M x M d ) and the Fourier multiplier (1 — x(D x )) A” 1 
is bounded on L°°(M d ), because its symbol cr(£) = —|£|~ 2 (1 ~~ x(0) is smooth in 
W l and has (inverse) Fourier transform in L 1 (M d ); indeed repeated integrations by 
parts give |a:| fc a(a;) G L°° for every k > d — 1. 

We now observe that the potential Vo(t, x) falls in the class considered in [11,12], 
and the corresponding propagator Uo(t, x) was proved in [12, Theorem 3.1] to satisfy 
the estimates 

(46) \\U 0 (t,s)\\ L 2^ L 2 <C 
for \t — s| < T. 

The perturbed original Cauchy problem can be written in integral form as 

u(t) = Uo(t, s)u 0 — ifT 1 J U 0 (t, t)Vi(t, • )u(t ) dr. 

This is a Volterra-type integral equation, which by (45) and (46) is easily seen to 
have a unique solution in C([s — T, s + T]; L 2 (M d )), for every T > 0. By (46) it is 
also clear that the propagator U(t,s) satisfies (42). ■ 

We are ready to prove our main result. 

Proof of Theorem 1.1. Consider the operator i? (0 ^(t,s) defined by 

R (fl) (t,s) = E<®(t,s)-U(t,s). 

By Proposition 4.1 and (39) we can write 

i? (0 ^(t, s)f = —ih~ l J U(t, r)G ( ' 0 ' ) (r, s)f dr. 

Using (41) and (42) we deduce that 

(47) \\R(°\t,s)\\ L ^ L i<C(t-s) 2 
for 0 < \t — s\ < 5. 

We then proceed as in [12, Lemma 3.2], where the analogous result for smooth 
potentials was proved. 
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Namely, we have to estimate the L 2 —y L 2 norm of the operator 

E^\n,t,s) - U(t, s) = £(°)(Ml-i)£ ( 0 ) (*L-1,*L- 2 ) • • .^ ( 0 ) (ti,s) 

= (U(t, t L . i) + R^°\t, t L - 0 ) .. . s) + R w (t 1, s)) - C/(f, s). 

We can expand the right-hand side, grouping together the consecutive factors “of 
type U ” (using the evolution property of the propagator), and then estimate the 
L 2 —> L 2 norm of each ordered product by (42) and (47). After that it is an 
algebraic matter to handle the sum which arises and express the result in the 
desired form. For full details we refer to the proof of [12, Lemma 3.2], which can 
be repeated essentially verbatim in our framework (the seminorms || • || m which 
appear there have to be replaced by the norm of bounded operators on L 2 (M d )). 
This gives the desired formula (6). ■ 


5. The case of higher order parametrices 


In this section we assume higher regularity on the potential and we prove a 
stronger convergence theorem. Namely, let N G N, N > 1. Suppose that 

V(t,x ) satisfies Assumption (B) in Section 3 with k = d + 1 + N([d/2] + 3). 


Therefore we will apply the results of Section 3 for this value of k. 

Let 5 be the constant appearing in Proposition 3.5, so that the function S(t, s, x, y ) 
is well defined for 0 < \t — s\ < 5. 

In order to construct higher order parametrices, consider the functions aj(t, s, x, y ), 
j = 1, 2,..., N, defined by the formulas (cf. [11, Formulas (3.4) and (3.5)]) 

(48) ai(t, s, x, y) = exp ( - ^ J (r - s)A x o;(t, s, x{r),y) dr) 
and 

(49) 


a j(t, s, x, y) = --ai(t,s,x,i/) 


A x aj-i(T, s, x(r),y) 
a i(r, s, x(r),y) 


dr , 


for j = 2where co(t,s,x,y ) is the function defined in (35) and x{r) = 
n{T,s,y,ri{t,s,x,y)). 

We need the following result on the regularity of the functions aj(t,s,x,y). 


Proposition 5.1. For 0 < \t — s\ < 6 we have 

(50) ai(t, s,x,y) = l-(t- s) 2 r(t, s, x, y) 

and 


(51) l/oi(f, s, x, y) = 1 + (t - s) 2 r'(t, s, x, y), 

where the functions r(t, s,x,y) and r'(t,s,x,y) belong to a bounded subset of 

j^d+3+(N—l)([d/2]+3) 
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Moreover for 2 < j < N we have 
(52) \\aj(t,S, *, •) || ^jd+3+(N— j)([d/2]+3) — C\t 3 | 

for some constant C > 0. 


Proof. Let us now prove (50). It follows from Proposition 3.6 (with k — d + 1 + 
N([d/2 ]+3)) that A x cu(r, s, x , y) belongs to a bounded subset of fld+i+ N ([d/2]+3) (j^m^ 
for 0 < \t—s\ < 5, so that by Proposition 2.3 (with the map ( y , () ha ( x(t , s, y , C), y) 
in place of g ) we deduce that the same is true for 

A x u(r,s,x(T,s,y,(),y) 

(the function x(r,s,y,() is defined in Proposition 3.5). 

By Proposition 2.1 we see that the functions 


A x uj(t, s, x(t, s, y, [(r - s)/(t - s)]C), y) 
for0< |r — s\ < \t — s\ <5 belong to a bounded subset of 

-£_j-d-\-l-\-N([d/2]-\-3) — ([d/2]-\-l) — 1)([g£/2]+3) 





and again by Proposition 2.3 (with the map ( x , y) ha (y, ({t, s, x, y)) in place of g ) 
the same holds for 


A x uj(t, s, x{t , s, y, [(r - s)/(t - s)]((t, s, x, y), y) = A x u(t, s, x(r),y). 

' 

=x(t) 


Hence we obtain the estimate 

[ (r — s)A x u>(r, s, x(r),y) 


H d+3+(N-l)([d/2]+3) 


< C\t — s| 2 


By expanding the exponential in (48) as a power series in the Banach algebra 
Hlf +(N - 1 A ( V 2 ]+ 3 ) (M 2fi ) we deduce (50). The same arguments apply to l/ai(t, s, x, y), 
which gives (51). 

Finally, by the same arguments as above and induction on j one sees easily 
that the functions %(£, s,x,y) for 2 < j < N verify the estimates in (52) (as a 
basic step for j = 2 one uses that A x ai(r, s, x, y) belongs to a bounded subset of 
-^fz +1+ ^ Ar_1)( ^ //2 ' +3) (® 2d )> therefore A x ai(r, s, x(t), y) belongs to a bounded subset 

Of ^d+l+(JV-l)([d/2]+3)—([d/2]+l) ( R 2d) __ ^<I+3+(-/V—2)([d/2]+3) | 


We now define the amplitude 

N 


a'--\t,s,x,y) = ^2(ih Y 3 aj(t, s, x,y), N = 1,2,... 


l=i 
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and the corresponding oscillatory integral operator 

E (N \t,s)f(x) = 1 [ e lh ~ ls{t ’ s ’ x ’ v) a( N \t,s,x,y)f(y)dy. 

(2m (t - s)h) d / 2 J Rd 

It follows from Propositions 3.6, 5.1 and 2.6 that E( N \t, s ) is bounded in L 2 (R d ) 
and satisfies the estimate 

(53) \\E^ N \t,s)\\ L ^ L 2 <C 

for some constant C > 0, for 0 < \t — s| <5, possibly for a smaller value of S. 
Moreover, by Propositions 5.1 and 2.6 we have 

|| E (N \t,s) - E w (t , s)^l, 2 ^l 2 < C(t — s) 2 

so that, by Proposition 4.1, for every / G L 2 (R d ) we have 

lim E {N) (t,s)f = f 

t—ys 

in L 2 (R d ). 

On the other hand we also have (cf. [11, Formulas (3.9) and (3.11)]) 

(ihd t + ^h 2 A -V(t,x))E {N) (t,s)f = G {N) (t,s)f 

with 

G {N) (t, s)f = - [ e lh ~ ls ^ s ’ x ’ y) A x a N (t, s, x, y)(t, s, x, y)f(y) dy. 

(2m(t - s)h) d / 2 J Rd 

By (50) (if N = 1) or (52) (if N > 2) we see that the amplitude A x ajy satisfies the 
estimate 

|| A x O,]y (t , S, •, •) || H d + 1 ( R 2 d'j A G 1 1 s\ 

Hence Proposition 2.6 gives 

(54) \\G^ N \t,s)\\ L 2^ L 2 < Ch N+1 \t - s\ N+1 . 

We then consider the remainder operator 

R {N) (t, s)f := E^ N \t, s)f - U(t, s)f = -Ik' 1 1 ' U(t, t)GW(t, s)f dr. 

Using (54) and (42) we deduce that 

(55) \\R {N) (t,s)\\ L 2^ L 2 < Ch N \t - s\ N+2 

for 0 < \t — s| < 8. 

Consider now the composition 

E^in, t, s ) = EW(t, t L ^)E^ N \t L ^tL- 2 ) ■ ■ ■ E^ N \t 1 , s), 

for any subdivision Q : s — to < ti <...< tL — t of the interval [s,t] such that 

uj(Q) < 6. 
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We have the following result. 

Theorem 5.2. Assume that V(t,x ) satisfies Assumption (B) in Section 3 with 
k — d + 1 + N([d/2\ + 3), for some N > 1. 

For even/ T > 0 there exists C = C(T ) > 0 swc/i that for 0 < t — s < T and any 
subdivision R of the interval [s, i] with to(f2) < 5, and 0 < h < 1, we have 

(56) \\E {N \n,t,s) -u(t,s)\\ L 2^ L 2 <ch N uj(n) N+1 (t-s). 

Proof. As in the case considered in Section 4 (which essentially corresponds to 
N — 0), the desired result follows exactly as in [12, Lemma 3.2]; namely one 
expands the right-hand side of 

F (JV) (L>, t, s)—U(t, s ) = (U(t, tL-J+RWit, t L ^)) ... (U{t u s)+R w (t 1 , s))-U(t, s), 

and estimates the L 2 —>• L 2 norm of each factor. To this end one uses (42) for the 
factors “of type U" and (55) for the factors “of type R\ We refer to [12, Lemma 
3.2] for full details. ■ 
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